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In this paper we show that many properties of the baby skyrmions, which have been determined 
numerically, can be understood in terms of an analytic approximation. In particular, we show 
that the approximation captures properties of the multiskyrmion solutions (derived numerically) 
such as their stability towards decay into various channels, and that it is more accurate for the 
"new baby Skyrme modeV which describes anisotropic physical systems in terms of multiskyrmion 
fields with axial symmetry. Some universal characteristics of configurations of this kind are 
demonstrated, which do not depend on their topological number. 

1 Introduction 

It is known that the two-dimensional 0(3) a- model |1|] possesses metastable states which 
when perturbed can shrink or spread out due to the conformal (scale) invariance of the 
model [0, 0, f|. This implies that the metastable states can be of any size and so a fourth 
order in derivatives term, the so-called Skyrme term, needs to be added to break the 
scale invariance of the model. However the resulting energy functional has no minima 
and a further extra term is needed to stabilize the size of the corresponding solitons, ie 
a term which contains no derivatives of the field, often called the potential (or mass) 
term. In this case the field can be thought of as the magnetization vector of a two- 
dimensional ferromagnetic substance 0, and the potential term describes the coupling 
of the magnetization vector to a constant external magnetic field. As the extra terms 
contribute to the masses of the solitons their dependence deviates from a simple law in 
which the skyrmion mass is proportional to the skyrmion (topological) number and the 
two-skyrmion configuration becomes stable showing that the model possesses bound states 
1- 

In this paper we demonstrate that the simple analytical method used for the de- 
scription of the three-dimensional Skyrme model presented in can be used also to study 
various properties of the low-energy states of the corresponding two-dimensional cr-model 
when the parameters which determine the contributions of the Skyrme and the poten- 
tial term are not large. More precisely, it was possible to describe analytically the basic 
properties of the three-dimensional skyrmions for large baryon numbers ||, and so it is 
worthwhile to derive such a description for the two-dimensional 0(3) cr-model as well. In 
general, such analytical discussions of soliton models are useful as they lead to a better un- 
derstanding of the soliton properties. The two-dimensional 0(3) cr-model is widely used to 
describe ferromagnetic systems, high-temperature superconductivity, etc and so the results 
obtained here can be useful for the understanding of these phenomena. 

Our method is based on the ansatz introduced in || and is accurate for the so- 
called "new baby Skyrme modeV which describes anisotropic physical systems. In fact, its 
accuracy increases as the skyrmion number n increases, and this method allows to predict 
some universal properties of the ring-like configurations for large n, independently on its 
particular value. Although such models are not integrable, in the case where n is large the 
u new baby Skyrme modeV appears to have the properties of an integrable system. 



2 Near the Nonlinear 0(3) cr-Model 



The Lagrangian density of the 0(3) a-model with the additional terms introduced and 
discussed in || [7], |J is:Q 

Here d a = d/dx a ; x a , a = 0,1,2, refer to both time and spatial components of (t,x,y); and 
the field n is a scalar field with three components n a , a = 1,2,3, satisfying the condition 
n 2 = nl+nl+Til = 1. The constants g, e are free parameters, ie g 2 has the dimension of energy. 
It is useful to think of g 2 and 1/ge as natural units of energy and length respectively. The 
first term in ([I]) is familiar from er-models; the second term, fourth order in derivatives, 
is the analogue of the Skyrme term; while the last term is the potential term. In fact, 
the potentials for the "old baby Skyrme model" (OBM) and the "new baby Skyrme modeV (NBM) 

describing anisotropic systems are given by 

Vqbm = (1 ~ n 3 ) , 

V mM = \i± 2 {l-n 2 ) (2) 

respectively, and ju has the dimension of energy, so defines a second length scale in our 
model. Evidently, V^ BM < V OBM at fixed value of /x. 

In three spatial dimensions the Skyrme term is necessary for the existence of soliton 
solutions, but the inclusion of a potential is optional from the mathematical point of view. 
Physically, however, a potential of a certain form is required to give the pions a mass 
0. By contrast, in two dimensions a potential term has to be included in the above 



Lagrangian in order soliton solutions to exist. As it has been shown in |T0[, the different 
potential terms give quite different properties to the multiskyrmion configurations when 
the skyrmion number is large. Our analytical treatment here supports this conclusion, as 
shown in sections 3-5. 

We are only interested in configurations with finite energy, so we define the configu- 
ration space to be the space of all maps ri-.R 2 —>S 2 which tend to the constant field (0, 0, 1) 
(so-called vacuum) at spatial infinity 

lim n{x) = (0,0,1). (3) 

| cc I — >oo 

Thus every configuration fi may be regarded as a representative of a homotopy class in 
^(S 2 ) = Z and has a corresponding integer degree of the form 

deg [n] = -!- / d 2 x e ab n {d a n x d b n) . (4) 

The vacuum field is invariant under the symmetry group G = E 2 ~x SO(2) iso x P, where E 2 is 
the Euclidean group of translations and rotations in two dimensions which acts on fields 
via pull-back. SO{2) lso is the subgroup of the three-dimensional rotation group acting on 
S 2 which leaves the vacuum fixed. [We call its elements iso-rotations to distinguish them 
from rotations in physical space] . Finally P is a combined reflection in both space and the 
target space S 2 . 

We are interested in stationary points of deg[n] ^ 0, and the maximal subgroups of 
G under which such fields can be invariant are labelled by a nonzero integer n and consist 

1 The first few paragraphs of this section follow quite closely to |3, H] and are included to make the paper more selfconsistent. 



of spatial rotations by some angle a e [0, 2n] and simultaneous iso-rotation by -na. Fields 
invariant under such a group are of the form 



n\ = sin /(f) cos(n^ 



n>2 = sin /(f) sin(n(j 



713 = cos /(f) 



(5) 



where (f , 4>) are polar coordinates and /(f) is the profile function. Such fields are the ana- 
logues and generalizations of the hedgehog fields in the Skyrme model. This parametrisa- 
tion which involves azimuthal symmetry of the fields assumes that all the skyrmions sit on 
top of each other while forming the multiskyrmion configuration. 
It is easy to show that the degree of the field (Q) is 



deg [ft] = n 



(6) 



ie equal to the azimuthal winding number n. 

The corresponding static energy functional connected with the Lagrangian (|1]) for 
the OBM and NBM is equal to 



E d{n) OBM 



n 2 sin 2 / 
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n 2 f' 2 sin/ 2 
V 2 
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(7) 
(8) 



respectively. In (^) and (|8|) the length ( v /ge/l)~ 1 has been absorbed so that the scale size 
of the localized structures is a function of the dimensionless spatial coordinate r = y/gejlr 
while the dimensionless parameter a = \ijge. becomes the only nontrivial parameter of the 
model. Finiteness of the energy functional requires that the profile function has to satisfy 
the following boundary conditions: /(0) = n and /(oo) = 0. 

By setting cj> = cos / into ([7p the energy functional becomes 
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(9) 



and a similar expression for E cl (n) NBM . Parametrizing the field 0, using the ansatz introduced 
in B for the description of the three-dimensional skyrmions, as 
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leads after integration with respect to r to the following analytic energy expressions 
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Here p and r„ are parameters which still need to be determined by minimizing the en- 
ergy. In fact, r n corresponds to the radius of the nth-soliton configuration. Remark: The 
following Euler-type integrals have been used for the derivation of flTT| ) and (|T2J), see also [BJ 
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It can be easily proved that the minimization of the energies ([11]) and (|T2| ) implies that 
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and so the minimum value of the energies is equal to 
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It is obvious that the energy contributions of the Skyrme and the potential term are 
equal due to (|I4D, which is in agreement with the result obtained from Derrick's theorem. 
Equations (|15|JT^) provide an upper bound for the energies of baby-skyrmions for any value 
of p. To get lower bound, we should minimize the right-hand sides of (|I5|JI^) with respect 
to the parameter p. In what follows we investigate various cases which correspond to 
different values of the only nontrivial parameter of the model, a. 

First consider the case where a <c 1, ie for very small values of the model parameter. 
Observe that for a = the ansatz (|TDp is a solution of the model for p = 2n, which implies 
that p — > 2n as a — > 0. In this case, due to (0), the radius of the multiskyrmion configuration 
increases with n: (r™ m )o BM ~ rt 3/2 and (r™™) 2 BM — n 2 . Moreover, the configuration consists 
of a ring of thickness: 5 — Ar n /p thus <5 bm ~ 2n _1/4 and 8 KBM ~ const. Remark The ring 
thickness is determined as the difference of the values of <j> inside (which is equal to -1) 
and outside (which is equal to +1) of the ring (ie d<f> = 2) divided by its derivative at r = r„ 
where, due to (|10|), <f>(r n ) = and so 0'(r„) = p/2r n . 

T2l as solutions of the 



This kind of magnetic solitons have been observed in [11 



Landau-Lifshitz equations defining the dynamics of ferromagnets. [Note that the static 
solutions of the baby Skyrme model and the Landau-Lifshitz equations are related.] In 
general, 4> given by (|10D for p = 2n is a low-energy approximation of multiskyrmion config- 
urations (for n > l), since for n = 1 the corresponding energies given by ( |15D and ([16]) are 
infinite. Indeed, it is a matter of simple algebra to show that 
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For large n, the energies take the asymptotic values 



E c i{n) c 



Anng^l 1 



T a|, 



E c i(n) 



Anng z 1 



(18) 



Note that the energy of the OBM per unit skyrmion number increases as n increases, 
while the energy of the NBM per skyrmion decreases with increasing n to become constant 
for n large. In fact, the energies given by (17) are the upper bounds of the multiskyrmion 
energies since the exact profile function corresponding to the minimum of the energy differs 
from that given by ([1C]) . 



3 Perturbation Theory for the Model Parameter 

In this section energy corrections up to second or higher order with respect to the model 
parameter a have been obtained. The corresponding energies for OMB and NBM can be 
written as 

E cl (n)=^g 2 [f{p) + ah{p)] (19) 

where f(p) and h(p) can be evaluated from ( p~5|) and (|T6|), respectively. By letting p = 2n + e 
and expanding the energies ( |T5| ) and ( pT6|) up to second order in e we get f(p) = n + e 2 /(8n) 
and h(p) = h + ehi where hi = (2n)' 1 j3h . In fact, the corresponding functions for the OBM 
and the NBM are given by 

■ 2 " 7 y/l-l/n 2 , /3oBM= ![ C ot( 7 r/n)-i + ^— , (20) 
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Minimization of (|19[) with respect to e implies that e mm = -Aanhx = -2af3h . At large 
values of n, the parameteres e and p = 2n + e take the values 

/2n , > [2n 2 /3- 1 
£ ( n )oB M - -anW y, e(") NB M - 2ay ^ (22) 

pWobm - 2n - an \lY> p(")nbm - 2n + 2a^^y — -. (23) 

For any a, as n increases, the effective power p(n) OBM becomes negative and the approach 
based on the assumption that e OBM is small is not self-consistent (also, see next section). 
On the contrary, for NBM, p(n) NBM ~ 2n as n increases which implies that our consideration 
is self-consistent in this case. In terms of ([19])-(^1]), the energy per skyrmion of the n-th 
skyrmion configuration takes the value 

Ed(n) h 2 h 2 (3 2 

- — = l + a a — r , 24) 
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which gives us 
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For large n, the energies fl24f) take the asymptotic values 
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Note that the energies of the two models behave differently when we consider terms of the 
second order in the model parameter, ie terms ~ a 2 . Indeed, for the OBM the contribution 
to the energy is linearly proportional to the skyrmion number n, while for the NBM the 
contribution decreases rapidly as the skyrmion number increases. This implies that the 
linear approximation in a is accurate for the NBM since the quadratic term becomes 
negligible for large n. Numerical results obtained for different values of a for the OBM and 
NBM are presented in Table 1 and Table 2, respectively. 

As we have pointed out earlier, our method cannot describe the one-skyrmion con- 
figuration since the corresponding energies become infinite. However, by setting p = 2 + s 
to ( |T5| ) and ([16]) and expanding all terms up to third order in e < 1 we obtain 
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where 7 has different value for each model, ie 
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Note, that when only terms up to second order in e have been considered, the corresponding 
energy (|27|) simplified to: E cl = Airg 2 + ^ + 2ay^) , and the minimum occurs at: e 1 = 
2(a/v / 3) 2 / 5 - Finally, the minimum of ( P?D occurs at 
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and corresponds to a shift of e 1 since higher order corrections of e have been considered in 
(P7|)- The energy of the one-skyrmion configuration is 
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Equation (|30| ) implies that, for a single skyrmion, the energy expansion in a is propor- 
tional to a power of a instead of being linearly proportional to a (which is the case for the 
multiskyrmion configurations with n > 2), while its convergence is worse than for multi- 
skyrmions, especially for the NBM. In fact, when a = 0.4213 the first two terms in (|30| ) are 
equal to 1.807, while the next order term lowers this value down to 1.44, which compared 
with the exact value 1.564 obtained from numerical simulations gives an error of 7%. Note, 
that our one-skyrmion parametrization gives the same energy for both models, when ex- 
pansions only up to the lowest order in a have been considered: the difference appears only 
in the term ~ a^^/e in ( |27|) . 
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a = 0.316 (num) 
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Table 1: Energy per unit skyrmion number (in Aitg 2 ) for different values of the parameter a for the 
OBM case where corrections of second order in a have been taken into account. The last two lines contain 
the exact results obtained from the numerical simulations of multiskyrmions (n > 2) with ring-like shapes 
and (n > 3) with shapes other than ring-like fiofl , respectively. For the first case, we have solved numerically 
the equations using the hedgehog ansatz (||). 
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Table 2: Energy per unit skyrmion number for different values of the parameter a for the NBM. 
The last line contains the exact results determined by the numerical simulations |jl^| of multiskyrmions with 
ring-like shapes when a — 0.4213, for n < 6 coinciding with ours. 



By looking at the results presented in Table 1 and Table 2 it is clear that our approxi- 
mate method gives the energy values which are quite close to the exact values obtained by 
numerical simulations, especially for the NBM. In particular, the difference between the 
exact and approximate energy for a = 0.4213 is less than 0.5% for n > 6. For smaller values 
of a the agreement between analytical and numerical results is even better. In evident 
agreement with (^), the energies of the NBM skyrmions given in Table 2 are smaller than 
those of the OBM skyrmions (see Table i) at the same values of the model parameters. 

Note that, for the OBM (when a is small) the energy per skyrmion of a multi- 
skyrmion configuration with n > 2 is smaller compared to the single skyrmion energy and 
therefore, these configurations are bound states, stable with respect to the decay into n 
individual skyrmions. On the contrary, the ring-like OBM multiskyrmions with even n 
(where n > 4) are unstable with respect to the decay into two-skyrmion configurations, 
while configurations with odd n (where n > 5) are unstable with respect to the breakup 
into a two- and a three-skyrmion configurations. In addition, Table 1 and (^0[) show that 
for any n (where n^l) there is an upper limit for the model parameter: a < a cr (n) above 
which the n-th ring-like skyrmion configuration can decay into n individual skyrmions. 

Let us consider the case of n = 3 in more detail. As it can be observed from the 
energies (|17[) and ( |30|) when a < 0.77 the ring-like three-skyrmion configuration is stable 
with respect to the decay into a single and a two-skyrmion configuartion since 

Ex +E 2 -E 3 ~ 1.611 a 4/5 -ott ( - 1 ) , (31) 



and this difference becomes positive when and only when 
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Corrections to the energy for the skyrmion configurations with n = 1, 2, 3 of the higher order 
in a, lead to smaller critical values a cr (n). 

Since our fields with axial symmetry (|5]) and ( pT0|) correspond to ring-like solutions 
of the Euler-Lagrange equations [[1] for a = 0, they have to be solutions of the corresponding 
equations also as a — > 0, ie when a takes values in a small region close to zero. [In fact, 
this region becomes more narrow as n increases as in this limit the expansion in a becomes 
less convergent]. On the other hand, the lattice-like configurations (tripole for n = 3, 
quadrupole for n = 4, etc.) are solutions of the equations when a > a cr (n) for given n 
|5|, [TU], |T3|]. However, the transition from the ring-like configuration to any other minimal 
energy configuration is a phenomenon which has not been studied in much detail yet and 
deserves further investigation. 

Finally, it should be stressed that, in contrast to the linear approximation, the 
quadratic approximation given by ( ^5|) does not provide an upper bound for the energy. 



4 Away from the Nonlinear 0(3) cr-Model 

In the general case, for arbitrary values of the parameter a and the skyrmion number n, 
soliton solutions can be obtained by minimizing numerically the energy ([15]) and (|T6| ) with 
respect to the variable p. This way an upper bound is obtained for the corresponding 
energies since the profile function is given by ([10]). 

For large a at fixed n (or for large n at fixed a), the expansion (^0|) is not self- 
consistent for the OBM. However, some analytical results can also be obtained in this case 
since (|I5|) for large a can be approximated by 
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Expansion, up to second order terms of (|33]) with respect to p, gives 



E cl (n) OBM ~ 4^ 5 2 ^ (l + p J (34) 

where c 2 = 2(7r 2 /3 - 1); while its minimization implies that p min ~ V3c^ = 3.71 and so the 
corresponding energy is 

Ed(n) OBM _ 4 . an ^ c ^ 1/4 
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Note that, by contrast with the results obtained near the nonlinear 0(3) a- model, for 
large a the parameter p is independent of the skyrmion number n. Also, for a > n the 
skyrmion radius is proportional to the square root of the skyrmion number: r„ ~ n 1 / 2 ; 
while the skyrmion thickness is given by: S ~ r n /p — n 1 / 2 , and so the ring-like structure of 
the configuration is not very pronounced. 



Direct numerical minimization of (^3|) with respect to p gives p min = 4.5 and the 
corresponding value of the energy is 

Ed(n) =155Qn (36) 

The energy, which has been obtained by solving numerically the Euler-Lagrange equation 
@, is: E cl (n)/A-Kg 2 = 1.333 an. The profile function corresponding to this solution is given 
by: cos / = ^.(r^ - r 2 ) + 2^- - 1 for r < r n and / = for r > r„. This solution is quite different 
from our parametrization fllPP and thus, the 16% difference between the exact and the 
approximate solution is understandable. 

To conclude, we recall that for NBM the parametrization ([10|) works well for arbi- 
trary large n, and its accuracy increases with increasing n, as illustrated in Table 2. 



Properties of the Skyrmions: Mean Square Radii, Energy Den- 
sity, Moment of Inertia 



Many properties of multiskyrmions can be determined using the ansatz (|10|). For example, 
the mean square radius of the n-th multiskyrmion configuration takes the simple form 



<r 2 > n = \ I drr 2 (f>' 
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= 2 v}, (37) 

where r n is given by ([Tip for the OBM and NBM. For small a, it was shown in section 2 
that p = 2n, which implies that the mean square radius becomes 



<r 2 > n <r 2 > n (38 ) 

OBM • / / \ /7> "nrm / / \ fn 



sin(7r/n)V3^' NBM nsin(7r/n)v / 3 



which is equal to n, 8^/3^3, ttVE, ■ ■ ■ and 87r/3, 2ny/E, ... for n = 2, 3, 4, . . ., respectively. 

For the NBM, even for a large enough value of the parameter a, the analytical 



formula ([14]) with the power p taken from (|23|) gives the values of < r 2 >„ NBM in a remarkably 



good agreement with those obtained in numerical calculations. E.g., for n = 3 the analytical 
result is \f< r 2 > 3 = 2.987, in natural units of the model l/(ge/j), to be compared with 2.872 
obtained numerically. This agreement improves with increasing n, and for n — 12 we have 
\/< r 2 >i2 ~ 10.92 to be compared with 10.85 determined numerically. A similar agreement 
between analytical and numerical results takes place for the mean square radius of the 
energy distribution of multiskyrmions (the 3D-case has been considered in detail in H). 

Note that the one-skyrmion configuration is (still) a singular case since ( |37l) is not 
defined for n = 1. However, as we have shown earlier, by expressing p = 2 + e and expanding 



14] ) in e we get r 2 =1 = ^2e/3 which leads to 



<r^>i=2i/- — - (39) 
V 3p m2n 

for e mm given by (|29|). So, our approximate method shows that as the model parameter 
tends to zero the mean square radius of the one skyrmion field tends to infinity since 
< r 2 >i~ a~ 1/5 ; while since < r 2 > NBM (n) = y/ri < r 2 > OBM (n) in this case the mean square 
radius is given by (p9|) for both models. 



The average energy density per unit surface element is defined as 

= St- ,40) 

with S ~ 2r n /n, see discussion after (|T6|). For the NBM, when n is large, ( fiO| ) takes the 
constant value 



^ NB M^w^Vi +a J (41) 

ie is independent of n. So, (^) represents the fundamental property of this kind of mul- 
tiskyrmions. On the contrary, for the OBM when ring-like configurations (which do not 
correspond to the minimum of the energy || [H|) are considered, the energy density in- 
creases with n like — yjn at small values of a. 

Another quantity of physical significance determining the quantum corrections to 
the energy of skyrmions is the moment of inertia which has been considered for two- 



dimensional models in [13|. In order to obtain the energy quantum correction of the 
soliton, due to its rotation around the axis perpendicular to the plane in which the soliton 
is located, we have to take the i-dependent ansatz of the form 

ni = sin /(f) cos[n(0 — uit)), = sin /(f) s'm[n(cf> — ujt)], 713 = cos /(f). (42) 

Then the w-dependence of the energy is given by the simple formula: 

E rot = (43) 
where 6j, the so-called moment of inertia, is given by |13[] 



0j(n) - g 2 n 2 J d 2 r sin 2 / (l + a/' 2 ) . (44) 



Using (|i0|) and the relations 



4 



[1 + {r/r n f] 



p 2 sin(27r/p) 



p > 2 



1 fn M*r [ x (r/r n ) 2p dr 



(1 - 2 ) 



16./ ' /■ ~ Jo [l + (r/r n ) p }\ 



6p 

we find that at large values of n the moment of inertia simplifies to 

2 / 2n anp \ 



i, p>0 (45) 



<-),(„) ^Wn^(^ + pj ,4(>) 



which holds for any multiskyrmion configuration described by ansatz (0), for both models. 
For small values of a, letting p = 2n and taking r 2 given by (|I^) we find that 



2 _X 



©./Hnbm - Wnrl ( 1 + a\/ § ) , (47) 



which implies that the moment of inertia, for large n, is 



Qj(n)~E cl (n)r 2 n , (48) 

in agreement with simple quasi-classical arguments for the thin massive ring. Similar 
quasi-classical formulae have been obtained for the three-dimensional skyrmions (see, for 
example, Ref. || [14]]) where the moment of inertia was shown to be given by 0j = 2M B r%/3 
for large baryon numbers; this expression is valid for a classical spherical bubble with the 
mass concentrated in its shell. 



6 Conclusions 

In this paper we have presented an analytical approach for deriving approximate expres- 
sions to skyrmion solutions of the two-dimensional 0(3) cr-model. These approximations 
are very accurate for small values of the parameter a which determines the weight of the 
Skyrme and the potential term in the Lagrangian. For other values of the model parame- 
ter we have performed some numerical calculations and then combined them with further 
analytical work to investigate the binding and other properties of multiskyrmion states. 

Two models have been studied: the "old baby Skyrme modef and the "new baby Skyrme 
modeV which differ from each other in the form of the potentials @. For both models the 
a dependence of the energy of a single skyrmion differs from the cases where topological 
number n > 2. For the OBM, when a is small, the n = 3 skyrmion configuration is stable 
with respect to the decay into a single skyrmion and a two-skyrmion configuration while 
the ring-like multiskyrmion configurations with n > 4 are unstable with respect to the 
breakup into two- and three- skyrmion configurations. 

For the NBM, on the other hand, the hedgehog multiskyrmion configurations con- 
sidered in [|10| and here describe bound states since the energy per skyrmion decreases as 



the skyrmion number increases. We note that the results obtained for the NBM are simi- 
lar to the ones obtained for the three-dimensional model studied in [pj. In both cases the 
energy per skyrmion decreases as the skyrmion number increases. The three-dimensional 



skyrmions obtained using the rational map ansatz [I5fl, for large n, have the form of a 
bubble with energy and baryon number concentrated in the shell. The thickness and the 
energy density of the shell (which is analogues to the thickness of the ring in the two- 
dimensional case) do not depend on the skyrmion number [BJ. Similarly, in this paper we 
have shown that the two-dimensional baby skyrmions of the NBM, for large n, correspond 
to ring-like configurations with constant thickness and constant energy density per unit 
surface of the ring. The building material for these objects is a band of matter with con- 
stant thickness and average energy density per unit surface. Thus the baby skyrmions can 
be obtained as dimensional reductions of the three-dimensional skyrmions at large n; while 
the three-dimensional skyrmions can be derived from the two-dimensional baby skyrmions 
as dimensional extensions. 

In || it has been concluded that the Casimir energy, or quantum loop corrections, 
can destroy the binding properties of the two-skyrmion bound states. It would be worth 
to investigate the validity of this argument for the two- and three-skyrmion bound states 
of the NBM. Another interesting question is to determine to what extent the region of 
small enough a is of importance from the point of view of physics. For large a the method 
overestimates the skyrmion masses for the OBM but is accurate for the NBM, especially 
for large n. 



The existence of bound states of the three-dimensional skyrmions has rich phe- 
nomenological consequences in elementary particles and nuclear physics. It suggests pos- 
sible existence of multibaryons with nontrivial flavor, strangeness, charm or beauty; more 
details are given in ]14[] and references therein. Similarly, the existence of bound states of 



two-dimensional baby skyrmions with universal properties in the NBM, which describes 
anisotropic systems, can also have some consequences for the condensed state physics, 
which would be worth to investigate in detail. 

VBK is indebted to G. Holzwarth for drawing his attention to the paper ||; his 
work is supported by the Russian Foundation for Basic Research, grant 01-02-16615. TI 
thanks the Nuffield Foundation for a newly appointed lecturer award. 
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